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Page 11, line 7, symbols should read: 

\ = q b 3 t) a q e 
Page li, line 8, symbols should read: 

M z = q b 3 | a q a 
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W s = q b 3 k dq 2 



NATIONAL ADVISORY C0MHITTE3 FOR AERONAUTICS. 

TECHNICAL MEMORANDUM NO. 488. 

DETERMINATION OF THE AIR FORCES AND MOMENTS 
PRODUCED BY THE AILERONS OF AN AIRPLANE.* 
By 0. Wieselsberger and T« Asano. 

I. 

The quantitative knowledge of the air forces and moments 
produced by the deflection of the ailerons of an airplane is 
of considerable importance to the airplane builder for maneu- 
verability investigations, as well as for the strength calcula- 
tion of the wings. I n a work previously published by one of 
the authors, full particulars relative to the method employed 
in computing these forces and moments were given, and the re- 
sults of numerical calculations obtained under certain assump- 
tions were also presented.** On account of the great number of 
variables affecting this problem, the results therein obtained 
are not applicable in all cases. On the contrary, they are 
valid only for a definite aspect ratio and a definite slope of 
the lift curve against the angle of attack; or, speaking more 
precisely, they are valid only for a certain nondimensional 
parameter p obtained by division of both these quantities. 

*"Bestimmimg der durch die Querruder eines Tragflugels erzeug- 
ten Luftkrafte und Moments," Seitschrift fur Flugtechnik und 
Motorluftschiffahrt, July 14, 1988. 

**C. Wieselsberger, "Theoretische Untersuchungen liber die Quer- 
ruderwirkung beim Tragflugel." Report of the^ Aeronautical Re- 
search Institute, Tokyo Imperial University, No. 30, 1937. 
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The exact determination of these aerodynamic quantities for 
other values of the parameter "p would require further quite 
extensive numerical computations, similar to those made in the 
above-mentioned work. It would therefore be desirable to de- 
rive approximate results for other values cf the parameter from 
the" above-mentioned results. A rough approximation method for 
another lav/ of lift has already been given in a note in the 
above-mentioned work. Subsequently it was shown that a consid- 
erably better method can be given, by means of which the mc.gni- 
tude of the air forces and moments may be determined with suffi- 
cient precision from the present results for all values of p 
occurring in practice. There is no object in carrying the ap- 
proximation too far, since, as a rule, the actual lift distri- 
bution at the wing tips shows deviations from the theoretical 
distribution, so that an approximate evaluation of the distri- 
bution may be regarded as satisfactory. After a few brief 
remarks on the fundamentals of the exact computation, the method 
will be so presented that the lift distribution for deflected 
ailerons may be determined for other values of the parameter 
p from the results already obtained. Coefficients will then 
be given in the form of diagrams and numerical tables, from 
which the desired forces and moments can be easily obtained 
by substitution in the given equations. 
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II. 

In regard to the details of the exact computation we must 
refer to the above-mentioned work. The method employed was 
essentially derived from S. Trefft^.* It is also found in H. 
Glauert's bock** in a somewhat different form. Here we can 
only mention that the lift distribution is expressed by a 
Fourier series, the coefficients of which are so computed from 
a system of linear equations as to satisfy the conditions of 
Prandtl l s wing theory lying at the basis of the whole computa- 
tion. The fundamental equation, which expresses the require- 
ments of the wing theory, must, strictly speaking, be satisfied 
at all points along the wing span. In the present case we have 
confined ourselves to 16 points. Due to certain characteris- 
tics of symmetry, it is only necessary to carry out the calcu- 
lation for one wing. The calculation then consists in the 
solution of eight linear equations with eight unknown quanti- 
ties. In this connection we must also notice that the discon- 
tinuous change of the angle of attack, as is actually the case 
with aileron deflections, is represented by a transition re- 
gion of finite length. Therefore we do not obtain an exact 
solution, but one which gives, nevertheless, sufficiently ac- 
curate results for practical purposes. The first task consists 

♦"Prandtlsche Tragflachen- und Propeller-theorie , ff Zeitschrift 
fur angewar.dte Mathematik und Mechanik, 1931, p. 206. 

**"The Elements of Aerofoil and Airscrew Theory," Cambridge, 
1938 * 
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in the evaluation of the lift distribution* From this, we can 
then determine without difficulty all that is of value regard- 
ing the forces and moments acting on the wingi 

III. 

We assume a rectangular wing in straight flight. It is 
understood that all angles of attack will be measured from the 
position at which the lift disappears, and it is assumed that 
the lift is a linear function of the angle of attack. The lat- 
ter statement is practically true, so long as the air flows 
smoothly along the wing profile, that is, up to the point of 
separation on the suction side. Our results, therefore, are 
valid only for angles of attack at which no separation of the 
flow takes place. The angle of attack of the middle section 
of the wing (where there is no aileron to alter its shape), is 
denoted by a. At the place where the aileron is deflected 
we shall refer the angle of attack to the original chord. In 
order to be able to retain the form of presentation of the curves 
in Figure 10, the lift coefficients, as plotted against the 
angle of attack, differ from those of the original work. For 
any given aileron deflection or settings, we consider the al- 
tered wing contour as a new wing profile, the aerodynamic char- 
acteristics of which we assume to be known, (Experimental re- 
sults of this kind are found, e.g., in the works of E. A. Lan- 
dells, W. G. Cowley, and E. A. Griffith, "Tests of New Forms of 
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Wing Sections with Flops Extending Over the Whole Length of the 
Aerofoil. 11 M.P.L. Report 1914/15; and of G. J, Higgins and 
E . H. Jacobs, "The Effect of a Flap and Ailerons on the N.A.C.A. 
M-5 Airfoil Section." N.A.C.A. Report No. 260, 1927. The actu- 
al aerodynamic angle of attack of the altered profile will "be 
indicated by a + a q or a - a q . If the altered wing is 
turned as a whole, with constant aileron setting, from its nor- 
mal position to the position at which the lift disappears, in 
both cases we shall consider the angle of attack of the un- 
broken profile as measured with respect to the chord or some 
other fixed line of reference. In contrast with the profile 
of the middle section, the broken profile thus has an aerody- 
namic ally effective angle of attack equal to Oq, which we 
shall chiefly use in what follows (See Section VI regarding 
the determination of Oq from the experimental results). In 
the matter of opposite aileron settings, we assume that they 
are equal in both directions and that the air forces produced 
by the equal angles bear the same relation to the total air 
forces. 

Aileron deflection causes a change in lift distribution. 
It is evident that the resulting distribution can be consid- 
ered as made up of two parts, namely, the lift of a wing with 
constant angle of attack a over the whole span, and the lift 
of a wing with an angle of attack 0° at the middle add an 
angle of attack + Oq or -oq at the ends corresponding to the 
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length of the ailerons. Since the lift distribution for the 
first part (constant angle of attack along the whole span) is 
known, we may confine ourselves to the calculation of the sec- 
ond part. In addition to the most important practical case in 
which the ailerons are set in opposite directions, the calcu- 
lations were also made for the case of aileron settings in the 
same direction. Naturally the length I of the ailerons (Fig, 
1), is the dominating factor for the forces and moments. The 
calculations were made, therefore, for ailerons of four dif- 
ferent lengths. Indicating the span by b, the lift distri- 
bution was determined for the following values of Sl/b: 
2 I 

-jj- = 0,234 0,5 0,658 1,0 for opposite settings; 

■— - = 0.234 0.5 0.658 1.0 for like settings. 

The intermediate- aileron lengths can be derived from the 
curves. The aileron chord does not enter, as such, into the 
computation, since the altered profile may "be regarded as a new 
wing section with known aerodynamic ' characteristics. Its ef- 
fect is expressed by the difference a q in the angle of attack, 
Furthermore, the magnitude of the previously mentioned parame- 
ter p is a measure of the lift distribution. If t denotes 
the wing chord, then p is expressed in the following form: 

v Cjt * 

o x being defined by the expression c a = 2 c x a' c being 
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the lift coefficient end a' being the effective angle of at- 
tack 'expressed in circular measure. c x is, therefore, nothing 
more than a measure of the slope of the lift curve against the 
angle of attack for infinite span, being identical with 



However, since model experiments are usually conducted with 
wings of finite span, the value of o x for infinite span must 
be derived therefrom. For the case in which the experiments 
were conducted with an elliptical wing contour, this conversion 
becomes simple. According to the wing theory, c x acquires 
the following value: 



F D stands for the wing area, b Q for the span, and a 
for the geometric angle of attack, corresponding to the lift 
coefficient c a , measured from the attitude of zero lift. 
This formula is not strictly correct, when a rectangular wing 
is used for the experiments. The error accepted in the compro- 
mise amounts, however, to only a few per cent, so that the for- 
mula is sufficiently exact for most practical purposes.* Theo- 
retical researches on the lift of wing surfaces gives for c x 

a value of about rr.* * Der iving c x by the a b ove f ormula ,_f gom_ 

*" Concerning the Exact Calculation of c x from the Experimen- 
tal Results of a Rectangular Wing," K. Glauert's book, p. 145 ff. 

**Von Mises, "Zur Theorie des Tragf lachenauf triebes, " Zeit- 
schrift fur Flugtechnik und Motorluf tschif f ahrt , 1917, p. 157. 




appropriately taken from experimental data. 



°i = 



3a- 
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experimental data, one obtains, for the usual profile forms, 

values between 2.5 and tt . The original calculations of the 

lift distribution were made for o x = tt and for an aspect 

ratio ^ = 2 tt , so that the parameter p took the value of 
t 

4. The lift distributions obtained for this case for the given 
aileron lengths are shown in Figures 2 and 3 for opposite and 
like aileron settings on both sides of the wing. These curves 
show only the additional lift produced by the aileron deflec- 
tion. The total lift distribution is found by adding this lift 
distribution to that of a wing with constant angle of attack 
over the whole span corresponding to the angle of attack of the 
middle wing section, which is already known. The lift distri- 
bution for a rectangular wing of aspect ratio 2 tt and constant 
angle of attack Oq over the entire span is indicated by dotted 
lines in both figures, 

IV, 

If we now undertake to determine the lift distribution for 
other values of the parameter p, we must, strictly speaking, 
again find the troublesome solution of the ^system of eight lin- 
ear equations, already mentioned. For another value of p the 
ordinates of lift distribution given in Figures 2 and 3 do not 
vary proportionally, but the shape of the whole lift curve var- 
ies; or, expressed mathematically, the different parameters 
p of corresponding distribution curves cannot be derived from 
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one another by a common transformation. We presume, however, 
that in case the value of the parameter for the calculation in 
question does not differ too much from 4 (which holds true for 
practical applications), tlje form of the distribution curve 
will not differ materially from the calculated curve. For oth- 
er values of p the ordinates of the curves would then he eas- 
ily increased or decreased in a certain ratio, i.e., the dis- 
tribution curves could be looked upon as related to one another. 
In order to prove the correctness of this assumption, the exact 

computation was carried through for the case of codirectional 

pi 

aileron deflections and for ~ = 0.5. Furthermore, the ex- 
act computation was carried out for a parameter p = 5 and 
compared with the corresponding results of the case p = 4. 
The result of this comparison is shown in Figure 4. Besides 
the two distribution curves, and B 5 ? of the lift on the 

semispan for p = 4 and p = 5 (the first being taken from 
Figure 3), there are also sketched the distributions A 4 end 
A 5 for constant angle of attack a q along the entire span for 
the same values of p. The ordinates of the curves B 4 were 
reduced in the ratio ^ of the ordinates of A 4 and A s . 
The distribution 3 5 * thus derived is represented by a dotted 
line. It can be seen that the latter does not quite coincide 
with the curve B s . We will, therefore, consider the distri- 
bution B B ! thus derived from B 4 as equivalent to the actual 
distribution B s and use it in the determination of the air 
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forces and moments. The resulting error is negligible for most. 

practical appl i cat ions . Since the ordinate ratio is some- 

h 4 

what different from place to place (the curves A 4 and A 5 them- 
selves being uncorrelated) , in the computation we shall at times 
use the ratio in the neighborhood of the maximum B 4 as indi- 
cated in Figure 4, 

V. 

From the lift distribution thus found, the additional 
forces and moments acting on the wing can now be determined. 
In the most important practical cases in which the aileron 
settings are in opposite directions, the following forces and 
moments act upon the wing: 

1. A rolling moment M x about the longitudinal axis, 

2. An additional induced dreg W x , 

3. A yawing moment M z about the vertical axis. 

The total lift in this case remains unchanged. The yawing 
moment M z is due to the fact that the unequal lift distribu- 
tion over the two wings produce, a different induced drag on 
each wing. The additional induged drag is caused by the con- 
siderable deviation from the most favorable lift distribution, 
which is known to be elliptical. 

Due to the complete symmetry, no yawing moment occurs for 
like settings of the ailerons. Pitching moments fall outside 
the field of our observations. Here we have: 
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4. An additional lift A x , 

5. An additional induced drag W 3 . 

As shown in the detailed study, these quantities may be 
expressed in the following manner; f 

For opposite aileron settings: 
Rolling moment M x = q b 3 £ aq 

Additional induced drag W x = q tP r\ aq 

Yawing moment M z = q b 3 C °q 

For like aileron settings: 

Additional lift A 1 = q b s \ cc q 

Additional induced drag W 2 - q b 3 k a q 

In these equations, besides the well-known quantities 
(q = dyriMiic pressure), there appear also the nondimensional 
coefficients £ , r\ , | 3 \, and k # These coefficients depend 
on the form of the lift distribution and are the result of in- 
tegrations over the span. The approximate values of these coef- 
ficients, deduced from the calculations for p = 4 with the 
help of the process given in Section IV, are reproduced in the 
curve groups of Figures 5 to 9. They are dependent on the ratio 
^± for different values of the parameter p; namely, for 
p = 4, 4.5, 5.0, 5.5, and 6. The coefficients £>••• * 
for all aileron lengths and for all values of p occurring in 
practice may be interpolated from these diagrams, and the de- 
sired moments and forces may be obtained in a simple manner by 
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substitution in the above formulas. The values of these coef 
ficients are also given in the following tables. 



1. Coefficients £ of the Rolling Moment 



21 

TT - ! 


0.234 


0.500 


0.658 


1.000 


p = 4.0 
4.5 
5.0 
5.5 
6.0 


0.047 
0.044 
0.041 
0.039 
0.038 


0.100 
0.092 
0.035 
0.079 
0.074 


0.114 
0.105 
0.096 
0.089 
0.084 


0.135 
0.125 
0.114 
0.105 
0.098 


2. Coefficients r\ of the Additional Induced Drag W a 


21 

b~ - 


0.234 


0.500 


0.653 


1.000 


p = 4.0 
4.5 

5.0. 
5.5 
6.0 


0.0462 
0.0402 
0.0353 
0.0318 
0.0297 


0.1114 
0.0939 
0.0806 
0.0696 
0.0607 


0.144 

0.121 

0.103 

0.0883 

0.0770 


0.198 
0.165 
0.140 
0.120 
0.104 


3. Coefficients i of the Yawing Moment M z 


21 _ 

"b 


0.234 


0.500 


•0.658 


1.000 

0.109 

0.0910 

0.0772 

0.0559 

0.0570 


p = 4,0 
4.5 
5,0 
5,5 
6,0 


0.0492 
0.0423 
0.0376 
0.0339 
0.0316 


0.0895 
0.0755 
0.0648 
0.0559 
0.0488 


0.09S5 
0.0613 
0.0690 
0.0592 
0.0520 
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4. Coefficients X of the Additional Lift A x 



21 
b " 


0.234 


0.500 


0.658 


0.826 


1.000 


p = 4.0 


0.130 


0.326 


0.442. 


0.587 


0.729 


4.5 


0.121 


0.299 


0.405' 


0.536 


0.663 


5.0 


0.113 


0.277 


0.374 


0.494 


0.607 


5.5 


0.106 


0.258 


0.346 


0.456 


0.559 


6.0 


0.100 


0.241 


0.324 


0.425 


0.519 



5* Coefficients k of the Additional Induced Drag % 



21 

"F - 


0.234 


0.500 


0.658 


0.826 


1.000 


p = 4.0 


0.056 


0.116 


0.137 


0.1^3 


0.178 


4.5 


0.048 


0.098 


0.115 


0.136 


0.147 


5.0 


0.042 


0.083 


0.098 


0.115 


0.123 


5.5 


0.037 


0,073 


0.084 


0.098 


0.105 


6.0 


0.033 


0.063 


0.074 


0.085 


0.090 


The quantities c l 


VI. 

and Oq 


for the 


determination of 



forces and moments for a given case will be taken from profile 
tests with full-length ailerons or flaps. The data obtained 
fey Higgins and Jacobs in the high-pressure tunnel, are here re- 
produced by way of illustration* The lift coefficients c n 
for different aileron settings (3, as well as the airfoil used 
(the N.A.C.A. M-6 with an aspect ratio of 6) are shown in Fig- 
ure 10* The aileron here extends along the whole length of the 
wing model, and the aileron chord is 30 per cent of the wing 
chord. It can be seen from the diagram that the lift curves 
run almost parallel to one another, that is, that the value 
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c i = o 5 „ 1S ver y nearly constant for various aileron set- 
d a a 

tings. From these curves we learn, furthermore, that the abso- 
lute vcJuc of the lift variation for positive and negative ail- 
eron settings, for the profile in question, is very nearly the 
same, so that the assumption made in this connection under Sec- 
tion, III holds good. However, it should be noted that a symmet- 
rical aileron effect cannot be expected for all profiles (i.e., 
equal lift variation on both sides for equal aileron settings 
in opposite directions), and that occasionally c x and a q 
may have different values for the two wings. In cases in 
which the asymmetry is not pronounced, fairly' good results are 
obtained by determining the air forces and moments for the two 
wings separately and simply adding the values thus obtained, or 
(which amounts to the same thing) by taking the arithmetical 
mean of the quantities c x and a q or of c/ and a q ' of 
the values found for both wings together. 

If, by way of example, we assume that the N.A.C.A. profile 
M-6 has an angle of attack of 3° to the direction of flight, 
as shown by the dotted reference line, the angle of attack 
measured from the attitude of zero lift would be about a « 4.5° 
If we further assume that the ailerons are deflected by 
0 = ±10°, we then see from the diagram that, for the angles 
of attack within the aileron working range, a + a q = +9.2° 
and a _ ctq = -0.2°. Hence the angle a q is 4.7°. From Figure 
10 we also learn that, for the customary form of lift-curve rep- 
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resentation (i.e., even though, for the deflected aileron, .the 
angle of attack is computed from the original reference line), 
the angle a q is only the angular difference between the inter 
sections on the angle- of -at tack axis of the curves (3 = 0° 
and p = 10°. 

Summary 

In an earlier work the lift distribution, produced by ail- 
eron deflections (in the same and in opposite directions) for 
a rectangular wing in straight flight, was determined on theo- 
retical principles for a certain case, and the air forces and 
moments produced by the ailerons were derived therefrom. In 
order to widen the field of ap'/xLicat ion, approximate values 
are here deduced from the former results, with the help of 
which the acting forces and moments may be calculated for all 
cases occurring in practice. The results are expressed by non- 
dimensional coefficients, which are presented in the form of 
diagrams* For any given case the corresponding coefficients 
are derived from these diagrams (by interpolation, when neces- 
sary), end, by inserting them in the given formulas, the de- 
sired forces and moments can be easily ascertained. 

Translation by 

National Advisory Committee 

for Aeronautics. 
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Fig. 2 Distribution of additional lift for different 
aileron lengths and opposite aileron settings. 



A 




Fig. 3 Distribution of additional lift for different 
aileron lenjfcF-e and like aileron settings. 




Fig. 4, 
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Figs. 5, 6, 7. 



0 .2 .4 .6 .8 1.0 

21/b 

Fig. 5 Coefficient £ of rolling moment. 



p=4.0 




0 



1.0 



.2 .4 .6 .8 

21/b 

Fig. 6 Coefficient r\ of additional induced drag for 
opposite aileron settings. 
,12r 




.2 .4 .6 .8 1.0 

21 /b- 

Fig.7 Coefficients i of yawing moment. 
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Figs. 8,9,10. 




p=4.0 
4.5 



Fig. 8 Coefficients 
6.0 * of additional lift 
for like aileron 
settings. 



.8 1.0 




1P=4.0 
4.5 



Fig. 9 



5,0 Coefficient k of 
5.5 additional induced 
6.0 drag for like aileron 
settings. 




Fig. 10 

Lift coefficient c a vs. 
angle of attack for 
different settings, P 
of full-length ailerons, 
on N.A.C.A.-M6 airfoil. 



16 24 " 32 40 48° 
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